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ABSTRACT 

We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, poly- 
tropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially 
in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of 
a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional insta- 
bilities, and is thus an unstable version of acoustic-surfacc-gravity waves. The Jeans mode dominates 
in weakly confined disks or disks with rigid boundaries. When the disk has free boundaries and is 
strongly pressure-confined, on the other hand, the mixed GI is dominated by the distortional mode 
that is surface-gravity waves driven unstable under own gravity and thus incompressible. We demon- 
strate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that 
inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective 
sound speed c e ff of acoustic-surface-gravity waves. We also find expressions for the gravity reduction 
factors relative to a razor-thin counterpart, appropriate for the Jeans and distortional modes. The 
usual razor-thin dispersion relation after correcting for c c ff and the reduction factors closely matches 
the numerical results obtained by solving a full set of linearized equations. The effective sound speed 
generalizes the Toomrc stability parameter of the Jeans mode to allow for the mixed GI of vertically- 
stratified, pressure-confined disks. 

Subject headings: Hydrodynamics — Instabilities — ISM: kinematics and dynamics — stars: formation 
— waves 
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1. INTRODUCTION 

Gravitational Instability (GI) plays a crucial role in 
structure formation in various astronomical situations 
ranging from growth of primordial de nsity fluctuations 
to star formation in disk galaxies (e.g.. lZerDovichlll97Ct 
iMcKee fc Ostrikerl [20071 ). GI in flattened systems is 
of particular importa nce in formation of giant clouds 
in disk ga l axies ( e.g.. IGoldreich fe Lvnden-Belll Il965at 
lElmegreeij 119871: iKim et all 120031) gas giant planets 
in pr otoplanetary disks (e.g.. iBosa 119971 : iDurisen et al.l 
|2007[ ). and bound clumps resulting from fragmenta- 
tion of shocked layers in star- forming regions (e.g., 
lElmegreen fc Ladalfl977l lElmegreenlll998[) . 

While real astrophysical disks are vertically stratified, 
most analytic studies often neglect the vertical degree 
of freedom and adopt an isothermal equation of state 
(EOS) for simplicity. In a rotating, infinitcsimally-thin 
disk with surface density Eo and sound speed c s , the lo- 
cal dispersion relation for axisymmetric waves with fre- 
quency lu and wavenumber k is given by 



27rGE fc 



(1) 



where kq is the epicycle frequency and G is the grav- 
itational constant (e.g., IGoldreich fc Tremainel 119781 : 
iBinnev fc Tremainel [20081) . If the disk is not rotating, 
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sound waves with k < kj = 2irGY,o/c^ would become 
gravitationally unstable. The Coriolis force arising from 
disk rotation stabilizes long-wavelength perturbations, 
making the disk Jeans-stable if the Toomre stability pa- 
rameter satisfies 



ttGSo 
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(2) 



for any k fe.g.. lToomrelll964l) . When Q T , C < 1, on the 
other hand, thermal pressure and rotation are unable to 
stop the collapse of over-dense regions. 

The thin-disk approximation would be valid as long 
as vertical motions are unimportant and the scale of in- 
terest is much longer than the disk scale height. Nev- 
ertheless, it has been well known that finite disk thick- 
ness makes some quantitative changes to the character- 
istic wa v elengt hs and critical Qt values. For instance, 
ILedouxl fl951) showed that waves in self-gravitating, 
non-rotating, isothermal disks with scale height Ho = 
Cs/(7r67S ) become unstable if kH$ < 1, which can be 
compared with the unstable co ndition kHg < 2 of the 
razor-thin counterpart (see alsolSimonlll965[) . For rotat- 
ing disks, IGoldreich fc Lvnden-Belll (|1965aD considered 
the effect of the finite disk thickness and found that the 
stability condition changes to Qt,c = 0.676. This de- 
crease of the critical Qt value is due to the dilution of 
self-gravity at the disk midplanc in a vertically stratified 
disk, which can approximately be treated by multiply- 
ing the gravity reduction factor T = 1/(1 + kHo) to the 
second term of the right-hand-side of equation (j]) (e.g., 
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lElmegreenlH987t iKim et alJl200l IKim fc OstrikCTl[2007h . 

While the results of GI in isothermal disks are infor- 
mative, there are several issues of GI that need clari- 
fication in more general situations. First of all, there 
are many classes of gaseous disks that do not have con- 
stant temperature in the vertical direction. Examples in- 
clude optically th ick regions of planet-forming protoplan- 
etary disks (e.g. iBell et al.l ll997t iD'Alessio etallll998t 
iBolev et"aLl 120071) and accret i on disks around compac t 
objects (e.g.. lLa Dousl 119941 : iLubow fe Qgilvid 119981 ). 
These disks have often been modeled using a polytropic 
EOS (e.g.. ILubow fc QgiMel 11991 INelson et al.l l2000l: 
IMamatsashvili fc Rice! I2010T ) . While behavior of var- 
ious waves in polytropic disks have been studied ex- 
tensively (e.g.. iLin et all 119901: iKorvcanskv fc Pringld 
fl995t ILubow fc Qgilvid 119981 : IQgilvie fc Lubowl 119991 ). 
GI of such disks has not bee n explored in detail. 
IGoldreich fc Lvnden-Belll (|1965aD calculated a stabil- 
ity criterion of a uniformly-rotating polytropic disk, 
but was limited to th e case with the adiabatic index 
7 = 2. iLarsonl ([1985D studied GI of polytropic disks 
with arbitrary 7 by using approximate scaling rela- 
tions instead of solving perturbation equatio ns accu- 
rately. Recently, IMamatsashvili fc Ricd (|2010( ) consid- 
ered vertically-stratified polytropic disks and argued that 
inertial modes rather than acoustic modes become un- 
stable to self-gravitating perturbations. In this work, 
we shall clarify using both analytic and numerical ap- 
proaches that it is the acoustic modes rather than the 
inertial modes that become unstable. 

Second, since the sound speed varies with height 
in polytropic disks, it is questionable what kind of 
average is the most appropriate for c s if one wants 
to use the stability condition ([2]). Several numeri- 
cal studies used the vertically-ave raged sound speed 
flMeua et alJl2005b IBolev et al.ll200(l or simply midplane 
value (jRice et al.ll2003ll2005l ) to evaluate stability of pro- 
toplanet ary disks, but the usag e of th ese has yet to be 
justified. Mamatsa shvili fc Ricd (|2010l ) attempted to re- 
solve this ambiguity by introducing Qm = fi 2 /(AnGpoo) 
as a three-dimensional stability parameter, where Q is 
the local orbital angular velocity and poo is the midplane 
density, but the connection between Qm and Qt is un- 
certain since the former does not depend on c s explicitly. 
In this work, we define the effective sound speed that 
well represents the modal behavior of acoustic waves in 
a thermally stratified disk and show that this naturally 
relates Qt to Qm- 

The third issue involves disk truncation by exter- 
nal pressure. Galactic dis ks are confined by ram 
pressure of infalling g as (e.g.. iDubois fc Tevssier] I2008T) 
or hot halo gas ( e.g., IGoldreich fc Lvnden-Belll Il965at 
iLee fc Hongj l2007f ). Thin shells produced by super- 
novae, stellar winds, or expan ding H II regions are 
usually bounded by shocks (e . g .. iDeharveng et all 120051: 
iChurchwell eTaTl 12001 120071 ). lElmegreen fc Elmegreenl 
(|1978h showed that pressure-confined disks become un- 
stable at scales smaller than the Jeans wavelength 
Aj = c 2 /(7rGSo) °f unbounded disks. The mass of 
fragments produced by such GI is less than the crit- 
ical Bonnor-Ebert mass, and thus may not necessar- 
ily experience gravitational r unaway. When the con - 
fining pressure is very strong, ILubow fc Pringlel ()1993l ) 



showed that the GI becomes essentially the same as 
that of an incompressible disk, which is in stark con- 
trast to the conventional Jeans ins tabili ty of compres- 
sional disks. lUmekawa et al.l (|1999() and IWunsch et al.l 
(|201CD confir med such prediction s by u sing numerical 
simulations. IBovd fc Whitworthl (|2005|) suggested the 
GI of strongly confined disks as a potential candidate 
for forming "free-floating" planetary-mass objects (hav- 
ing mass as l ow as ~ 0.003M^) detected in young star 
clusters (e.g.. IZapatero Osorio et aI]|2002D . More elabo- 
rate theoretical and numerical models inclu ding the ef- 
fects of shell expansion and deceleration (e .g., lElmegreenl 
Il989t llwasaki et al.ll2011t iDale et al.ll2011[). the rmal and 



chemical processes dHosokawa fc Inut suka 2006), b ound 
ary conditions (e.g.. lElmegreenlll989l: lUsami et al.l|1995h 
IDale et al.ll2009l ). magnetic fields (|Nagai et al.lll998[ ). etc. 



examined consequences of GI occurring under high exter- 
nal pressure. 

Despite these efforts, however, the physical mechanism 
behind the GI of strongly pressure-c onfined disks still 
remai ns controversial. For example, ILubow fc Pringlel 
(|1993f ) argued that the instability is due to a neutral 
mode that exists because external pr essure can h o ld the 
laye r to any distor t ed sh ape, while lElmegreenl (|1989l ) 
and lUmekawa et"aLl (I1999D claimed that the distortion of 
the surfaces exerts a pinching force that causes the dis- 
torted disk to collapse. On th e other hand, IWunsch et al.l 
(|2010t ) and IDale et al.l (| 2 1 If) considered the collapse of 
an oblate spheroid with unif orm density pressurized b y 
an ambient medium (see also lBovd fc Whitworthll2005[ ). 
and termed the enhanced GI "pressure-assisted grav- 
itational instability", although the physical processes 
responsible fo r it w ere not identified. Very recently, 
llwasaki et all ()2011l ) studied GI of pressure-confined 
shells around expanding H II regions taking into account 
boundary effects. They suggested from the vertical be- 
haviors of eigenfunctions that the restoring force of the 
unstable mode comes from surface-gravity waves. We use 
both numerical and analytic approaches to confirm that 
it is indeed surface-gravity waves that become unstable 
in strongly confined disks. 

In this paper, we investigate GI of rotating, vertically- 
stratified, pressure-confined ga s disks. We adopt simple 
polytro pic, rotating models of IGoldreich fc Lvnden-Belll 
(|1965al) and extend them to pressure-confined situations. 
This work also extends the case of non-rotatin g, isother- 
mal models of lElmegreen fc Elmegreenl (|1978l ) to rotat- 
ing, polytropic disks. Our objectives are four-fold. First, 
we wish to distinguish conventional Jeans instability 
from the GI arising from the surface distortion of strongly 
pressure-confined disks, and give a clear physical expla- 
nation for the latter that has previously been confusing. 
Second, we provide analytic expressions for the effective 
sound speed c e ff that best represents the modal behavior 
of sound waves in polytropic disks. Third, we find the 
gravity reduction factors caused by finite disk thickness 
for the Jeans and distortional modes of GI. We show that 
the numerical dispersion relations found by solving a full 
set of perturbation equations are well matched by equa- 
tion (fT]) provided that c s is changed to c e g and that the 
appropriate reduction factor is considered in the gravity 
term. Finally, we generalize the Toomre stability param- 
eter of razor-thin disks to vertically-stratified, pressure- 
bounded disks. While our disk models arc simple with a 
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polytropic EOS and without considering external gravity, 
some qualitative properties of GI clarified in this paper 
can hold true in more general situations such as disks 
with a barotropic EOS and/or with external gravity. 

The remainder of this paper is organized as follows. 
In Section [21 we introduce basic equations, construct an 
equilibrium model, and present linearized perturbation 
equations. In Section [3j we derive expressions for the ef- 
fective sound speed c c g that describe the phase speeds of 
fundamental modes in non-self-gravitating, non-rotating 
disks. In Section |4l we present the gravity reduction 
factors J 7 , and show that an approximate dispersion re- 
lation using c e ff and T is in excellent agreement with 
the numerical results obtained by solving the perturba- 
tion equations. For rotating disks, we define a stability 
parameter using c c e that is applicable also to the distor- 
tional mode. In Section[5l we summarize and discuss our 
results. 

2. FORMULATIONS 

We consider a rotating, self-gravitating, gaseous disk 
in vertical hydrostatic equilibrium and analyze its sta- 
bility to small-amplitude, axisymmetric perturbations. 
The disk is infinitely extended in the horizontal direc- 
tion, but truncated at a finite height by a tenuous ex- 
ternal medium with pressure p cx t- For local modes in 
the horizontal direction, it is advantageous to consider 
a local Cartesian reference frame whose center lies at a 
radius Rq and orbits the disk center with a fixed an- 
gular velocity fio = Q(Ro)- In this local frame, radial, 
azimuthal, and vertical coordinates are represented by 
x = R — Rq, y = Ro(4> — Clot), and z, respectively, 
and t he terms arising from curvature effect are ignored 
Ce.g iGoldreich fc Lvnden-Belll [19651 Uulian fc Toomrd 
ll966tlKim et al.ll2002[ ). Assuming that the shear rate 
q = — d In fl/d In R\n is uniform, the equilibrium back- 
ground velocity field is given by Uo = — qfloxe y and the 
local epicyclic frequency is Kq = (4 — 2g)f2§. The equa- 
tions of ideal hydrodynamics expanded in the local frame 
are 



0. 



(3) 



d\x 1 

— + u • Vu = — Vp - VV> + 2q£llxe x - 2fl x u, (4) 



V 2 ^ = 47rGp, 



(5) 



where p, u, p, and ip denote mass density, velocity, ther- 
mal pressure, and self-gravitational potential of the gas, 
respectively. To close the set of equations ([31)-©, we 
assume that the gas follows a polytropic EOS 



p = Kp 1 



(6) 



where K is an arbitrary constant and 7 is the polytropic 
exponent. If 7 is different from an adiabatic exponent 
of disturbances, the system under consideration exhibits 
complicated modal behaviors associated with co nvective 
motions. However, Mamats ashvili fc R ice (20lIJ showed 
that the convective modes do not affect the properties of 
GI much. In this work, we thus simply take 7 equal to 
the adiabatic exponent. 



2.1. Initial Equilibria 

We first explore equilibrium profiles of density and 
temperature in purely self-gravitating, polytropic disks 
confined by an external medium. The pressure con- 
finement model we adopt is a s traightforward extension 
of isot hermal disks studied by lElmegreen fc Elmegrecn 
(|1978| ) to polytropic cases. For hydrostatic equilibrium 
in the vertical direction, equations ([4]) and ([5j) require 



1 dp 
Po dz 



dipp 
dz 



—^0 = 47rGp , 



(7) 



(8) 



where the s ubscript "0" denotes the unpert urbed state. 

Following lElmegreen fc Elmegreenl (|1978t) . it is conve- 
nient to convert the vertical coordinate z to the dimcn- 
sionlcss variable p defined by 



1 



dipo 



AttGpqqHq dz ' 
where poo = Po(0) and Hq is the disk scale height 

K 



(9) 



(10) 



In Appendix [S] we show that equations ([7]) and ([5]) sub- 
ject to the symmetry condition dpo/dz — at z = yield 
the solution 

Po = Poo(l-M 2 ) 1/7 , (11) 

(see also lHarrison fc Lake! [1972). The corresponding lo- 
cal speed of sound is 



^lPo/po = c 2 s0 (l-p 2 ) 1 - 1 ^, 



(12) 



where c s0 = (7Poo/poo) 1/2 = (-fKPoo *) 1/2 is tnc mid- 
plane sound speed. 

Figure [T] plots equilibrium structures of self-gravitating 
polytropic disks with differing 7. For a softer EOS 
(smaller 7), mass density falls off more rapidly from the 
midplane to provide the pressure support against the 
gravity ~ -iirGpooz at \z\/Hq < 1. This in turn makes 
the gravity smaller at high-altitude regions, causing the 
density to decrease slower with height at \z\/Hq > 1. As 
7 increases, therefore, the density becomes more concen- 
trated toward the midplane. For 7 > 1, equation (|A9j) 
implies that the density becomes automatically zero at 
finite height 



kcutl = r(i-i/ 7 ) 

H 2 r(3/2-l/ 7 ) 



(13) 



where T is the Gamma function, indicating that disks 
with 7 > 1 arc truncated even without external pres- 
sure. The sound speed in a disk with 7 > 1 decreases 
monotonically with increasing z. Note that |z cut | = 00 
for 7 < 1. 

Now we allow for the truncation of a disk at height 
z = ±a due to an external pressure p oxt . The condition 
of pressure equilibrium at the boundaries requires 



A = ji{z = ±a) 



2 Pc 



ttGE 2 , 



-1/2 



(14) 
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2.5 



z/H 



Figure 1. Distributions of density po (top), local sound speed c s 
(middle), and vertical gravity go = dtpo/dz (bottom) in hydrostatic 
equilibrium of self-gravitating polytropic disks with exponent 7. A 
stiffcr equation of state (larger 7) leads to a flatter density profile 
near the midplane. Even without externa l pr essure, disks with 
7 > 1 arc truncated at a finite height (eq. 1131 V For 7 < 1, c s is 
an increasing function of z. 



(e.g. , lElmegreen fc Elmegreenl Il978f h Using equation 
(|A8[) . the column density of the bounded layer is given 

by 



S 



po dz = 2p 00 AH = 2p 00 H . 



(15) 



where H = AH$ is the effective half-thickness of the 
layer. The relative density drop at the boundaries is 
Poa/poo = (1 - A 2 ) 1 ^ where p 0a = p (z = a). Note that 
A = 1 when a = |-z cu t|, corresponding to no pressure 
truncation. As p ex t increases, the vertical sound cross- 
ing time (~ H/c s0 ) becomes smaller than the internal 
free-fall time (~ {2irGp 00 y 1/2 ). In the limit of A -> 0, 
the disk takes uniform density and H —¥ a regardless of 
7: in this case, self-gravity is so weak that the initial 
equilibrium is maintained by balance between p cxt and 
Poo- 

2.2. Perturbation Equations 

We apply small-amplitude perturbations to the initial 
equilibrium configurations described above. Equations 
©-([S]) are linearized to 



^+V-(p u 1 ) = 0, 



(16) 



^7- = ^rVpo - — Vpi - VV»i + qflou lx e y - 2O x u x , 
dt p 2 po 

(17) 
(18) 



V 2 i/>i = AirGp! , 



where the subscript "1" indicates the perturbed quanti- 
ties. 



We further assume that perturbations are axisymmct- 
ric and of plane- wave type 



qi(x, z, t) = q\{z) e~xp(iujt + ik x x) , 



(19) 



where q\ refers to any perturbed variable with frequency 
uj and radial wavenumber k x . Then, equations (|16 [) — (|18 [) 
can be reduced to four first-order ordinary differential 
equations 



#1 



dhi 2 

-r~ = w si* 

dz 

#1 , / 
dz 

dip'i ,2 , 
-dz~- k ^- 



+ipi)- — {h x - goi lz ) , 

I c s 

Ip'l , 



4nGp 



(20) 

(21) 
(22) 
(23) 



where £i z is the vertical Lagrangian displacement defined 
through u\ z = d^iz/dt, hi = pi/po is the perturbed 
enthalpy, and ip[ = dipi/dz. We take hi, ipi, and ip[ 
as four independent variables. 

The set of the perturbation equations (|20 ]) -([23 |) are 
to be integrated numerically along the vertical direction 
to yield a dispersion relation u> = u(k x ) that satisfies 
the boundary conditions (BCs) as well as the constraints 
of even-symmetry modes. To separate the distortional 
modes of GI from the conventional Jeans modes, we con- 
sider two distinct conditions at the boundaries: the rigid 
boundary that allows only the Jeans modes and the free 
boundary in which Jeans and distortional modes coexist. 
Appendix [B] describes the BCs we adopt. Appendix [C] 
presents the numerical method we follow to obtain full 
dispersion relations, which will be compared with ap- 
proximate dispersion relations in Sections [3] and 21 

For future purposes, we note that the total perturbed 
density p\ is a superposition of the perturbed density pn 
inside the disk due to wave motions and the perturbed 
density pi. s at z = ±0 due to the surface distortion: 

(24) 



Pi = Pi, 



Pi, 



where 



Pi. 



Poa£iz\8(\z\ - a), 



(25) 



for even-symmetry modes. Here, S is the Dirac delta 
function. The corresponding perturbations in the surface 
density can be written as Si = £1,; + Si, s . We define 



w 



Si, 



Si.j + Si. 



(26) 



as the fraction of the perturbed surface density inside 
the disk relative to the total perturbed surface density. 
We will show that w defined in equation (j2"f5]) is a key 
parameter that controls the relative importance of the 
Jeans modes to the distortional modes. 

2.3. Classification of Local Modes 

Ignoring waves arising from surface distortion, the per- 
turbation equations ([2!I]) -(f2"3" f give rise to two principal 
types of waves: inertial waves and acoustic waves, both 
modified by self-gravity The former is characterized by 
epicycle mot ions, while the latter is based on gas com- 
pressibility. iMamatsashvili fc Ricd (|2010l ) argued that 
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inertial modes, rather than acoustic modes, are strongly 
influenced by self-gravity to become unstable. In Ap- 
pendix [D] we classify two types of waves using the local 
dispersion relations and compare them with the numeri- 
cal results. We directly demonstrate that while the iner- 
tial modes are not much affected by self-gravity and thus 
remain stable, it is the fundamental acoustic modes that 
can become unstable in the presence of self-gravity. This 
makes sense since both acoustic waves and self-gravity 
rely on density perturbations, while inertial waves are 
incompressible in nature (e.g.. lLatter &: B albus 2009). 

Before exploring numerical dispersion relations, some 
physical insight on the fundamental acoustic waves in 
disks can be gleaned by integrating equation ([16)) and 
the horizontal component of equation (|17l) multiplied by 
po over z. With the help of equation (p~9|) . we obtain 

2 ,2 V + 2 / 07 N 

w = fc^o— + k , (27) 

where the hat indicates the density-weighted vertical av- 
erages (i.e., qi = J po qi dz /T,q) . Because hi, ^i, an d Si 
are interrelated with each other and generally depend on 
k x and u, one needs to solve equations (|2"0")) -(f2"3" |) numer- 
ically to obtain full dispersion relations. Nevertheless, 
equation (f2"T)l implies that lu 2 consists of three terms, 
each responsible for the effect of pressure (including ther- 
mal pressure and surface distortion), self-gravity, and the 
Coriolis force, as in equation (TTJ). This opens a possibility 

that by finding solutions of hi and ipi independent of lj 
in some suitable limits, one can obtain an approximate 
dispersion relation that matches the numerical results 
closely. We will focus on this in the remainder of this 
paper. 

3. NON-SELF-GRAVITATING FUNDAMENTAL MODES 

As mentioned before, the sound speed of polytropic 
disks varies with the vertical height, so that it is inter- 
esting to find an effective sound speed that best repre- 
sents modal behavior of waves propagating throughout 
the disk. Also of interest is the effect of surface distortion 
on sound waves in the case of the free boundaries. To 
address these issues, we in this section limit ourselves to 
the fundamental acoustic modes in non-self-gravitating 
disks. For simplicity, we ignore the effect of disk rota- 
tion, which merely adds Kq to u> 2 (see eq. [2"7]). 

3.1. Rigid Boundary 

In non-sclf-gravitating, non-rotating disks, equation 
(|27|) is reduced to 

- 2 = ^ = i^^kl = c 2 cS kl (28) 

where c g is the effective sound speed. When the bound- 
aries are rigid or A — 1, £i. s = and 

2 _ Jpic 2 s dz CgPohi dz 
J Pi dz J_ a p hi/cidz 

Equation (|29[) shows that excited waves are pure acoustic 
modes (p modes) whose restoring force is the thermal 
pressure. Obviously, c e g = c s $ in an isothermal disk. 
Note that we did not make any approximation in deriving 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 , 1 r 




j i i i I i i i i I i i i i I i , i L 



0.0 0.5 1.0 1.5 2.0 

k x H 

Figure 2. Dispersion relations of the fundamental acoustic modes 
in non-self-gravitating, non-rotating, vertically-stratified, poly- 
tropic disks with A = 1 and 7 = 0.8, 1, 1.5, 2. The frequency u> 
in the ordinate is norm alize d by c s nk x . Horizontal dotted lines 
mark c;?ff/ c s0 (= a) (eq. 1311 ) in the long wavelength limit. 

equations (f2"8")l and (f2T)|) except for ignoring the effects of 
self-gravity and rotation. 

Now consider the fundamental mode in the long- 
wavelength limit. Appendix [E] shows that the funda- 
mental acoustic waves have a simple solution h\(z) = 
constant in the limit of k x — > 0. Equation (|29|) is then 
simplified to 

f a On dz 
c 2 cii = c s 2 = t aF , for rigid BC, (30) 

J-a Po/ C t dz 

indicating that the effective sound speed corresponds to 
the density-weighted, harmonic mean of the local sound 
speeds. For polytropic disks, the ratio of the effective 
sound speed to the midplane sound speed is given ana- 
lytically by 

a = r~ 2 /4, = 2Fi(l |; A 2 f\ (31) 

where 2F1 denotes the Gaussian hypergeometric func- 
tion. Disks with 7 > 1 have a < 1 since the temperature 
decreases monotonically with height. When there is no 
external confining medium (^4 = 1), a varies smoothly 
from 3/2to2/7rfor0.5 < 7 < 2. Larger external pressure 
makes a closer to unity. 

To check the applicability of equation (|30|) for waves 
with finite wavelengths, we obtain numerical dispersion 
relations of pure acoustic waves in non-self-gravitating, 
non-rotating disks for A = 1 and 7 = 0.8, 1, 1.2, and 2, 
by following the procedures given in Appendix [Cj Fig- 
ure [2] plots the resulting numerical dispersion relations 
using solid lines, which can be compared with the ap- 
proximate dispersion relation (|28|) together with equation 
(|30j) shown as dotted lines. When 7 = 1, the dispersion 
relation is simply w 2 = c 2 Q k 2 . When 7 7^ 1, the numer- 
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Figure 3. Dispersion relations of the fundamental acoustic- 
surfacc-gravity modes in non-self-gravitating, non-rotating, poly- 
tropic disks with 7 = 0.8,1,1.5,2. The case with mild pressure- 
confinement with A = 0.8 (pext /"'GSg = 0.28) is shown. The 
frequency u is normalized by c a ok x . Solid lines plot the full nu- 
merical r esu lts, while dotted lines draw the approximate dispersion 
relation 1331 . 

ical dispersion relations progressively deviates from the 
analytic results with increasing k x H . When k x H < 1, 
however, the two agree within ~ 4% for 0.8 < 7 < 2. 
The agreement becomes better for disks with smaller A. 
Since GI occurs for waves with k x H < 1 when self-gravity 
is included, this suggests that c e g given in equation (|30|) 
is a good representative of the averaged sound speed. 

3.2. Free Boundary 
In the case of the free BCs, £i,«/£i » = 
P0a£iz\z=a/ Jq Pi,idz. Then, the effective sound speed 
in equation (p8|) becomes 



' S- a Pohi/c % s dz POa h la 



L cff 



/" a p hx dz g 0a J a p hi dz 



(32) 



where equation (|B3|) is used. Equation (|32[) indicates 
that waves in a pressure-confined disk with free bound- 
aries make use of two restoring forces: thermal pres- 
sure and vertical gravity. The former drives acous- 
tic waves, while the latter is responsible for surface- 
gravity waves. We term these mixed waves acoustic- 
surface-gravity waves. The fact that the effective sound 
speed has two terms, each arising from sound waves and 
surface-gravity waves, is completely analogous to longi- 
tudinal waves propagating in a distensible tube where 
the effective compressibility of a gas is the sum of the 
true gas compressibi lity and the distensibility of the tube 
(e.g.. lLighthil][l97cl . 

Since p Qa = p 00 (l - A 2 ) 1 /'* and g 0a = AwGp 00 AH , 
it is apparent that the acoustic term in equation (|32p 
dominates for A close to unity. On the other hand, the 
surface-gravity term becomes important in disks with 



Figure 4. Ratio of the perturbed surface density inside the disk 
(£1 i) to the total (£1 j + £1 s ) for acoustic-surface-gravity waves 
in the long wavelength limit as a function of A for differing 7. 

A <C 1. Appendix [E] shows that in highly confined 
disks, even fundamental modes have a particular solu- 
tion hi ex cosh(k x z). Equation (|3"2")l is then reduced to 



-eff 



1 



P0a/P 



Oil 



k r H 



= WC S 



for free BC, 



c s 2 g 0a H tanh(k x H) / 

with w given by equation (|26[) . The sonic contribution 
c s is valid for k x H < 1, as in the rigid BC case. The 
contribution from surface-gravity waves accounts for the 
reduced density at the surfaces. In the limit of A <C 
1, equations (|28|) and ([33]) recover the usual dispersion 
relation of surface-gravity waves 



= go a k x tanh(k x H), 



(34) 
(e.g., 



in an incompressible medium 

IGoldreich fc Lvnden-Be1illl965aD . 

Although equation ([3"51 is derived for a disk with cither 
A <C 1 or A ~ 1, we find that it is applicable for arbitrary 
A as long as k x H < 1. Figure [3] compares equation (|33[) 
(dotted lines) with the direct numerical solutions (solid 
lines) found by solving the full perturbation equations 
(with self-gravity and rotation neglected) for moderately 
confined disks with A = 0.8 and 7 = 0.8, 1, 1.5, 2. The 
discrepancies between the two are less than 0.5% for 
k x H < 1. This proves that c e ff given in equation (|3"3"|) is 
excellent in describing the phase velocity of the acoustic- 
surface-gravity waves in the long-wavelength regime. 

From equation (|33[) . it is now clear that the variable 
w defined by equation (|26l) is a weight factor repre- 
senting the relative importance of acoustic to surface- 
gravity terms in the dispersion relation of acoustic- 
surfacc-gravity waves. In the long-wavelength limit, 
equation (|32|) leads to 



w Q = w\ km ^ Q 



1 + "/a 



[l-A 2 fh 



2A 2 



(35) 
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When 7 = 1, w = 2A 2 /(1 + A 2 ) = (1 + p^/irGZZ)- 1 . 
Figure [4] plots wq as functions of A for disks with 7 = 
0.8, 1, 1.5, 2. As expected, surface-gravity modes dom- 
inate for smaller A and larger 7. From equations (|26|) 
and (|35[) . one finds that the mean amplitude of density 
perturbation inside the disk is given by 



So 



2A 2 a 



7a 

12U/7X 



(36) 



while Si jS /E = (1 - A z ) l n £,iz\ z=a l H at the disk sur- 
faces. Therefore, |Si,i/Si iS | < 1 for small A. This 
clearly demonstrates that the fundamental mode and 
its unstable version (upon inclusion of self-gravity) are 
incompressible in highly pres sure-confined disks (e.g. , 
Elmegreen fc Elmegreenl 119781 iLubow fe Pringli 119931: 



Nagai et al.lll998tlUmekawa et al.lll999[) . 

4. GRAVITATIONAL INSTABILITY 

In the preceding section, we derive expressions for 
the effective sound speed for fundamental acoustic 
or acoustic-surface-gravity waves in non-self-gravitating 
and non-rotating disks. In this section, we derive the 
gravity reduction factors caused by finite disk thickness 
that explain the numerical dispersion relations of GI 
fairly well. We begin by considering non-rotating, zero- 
temperature disks, and then include the effects of ther- 
mal pressure and rotation. 

4.1. Pressureless Disks 

Finite disk thickness is known to stabilize the conven- 
tional Jeans instability by diluting self-gravity at the 
disk midplane. As mentioned in Introduction, J- = 
(l + kxHo)^ 1 is often used as the gravity reduction factor 
for a n unbounded, exponenti a l disk with scale height H 
(e.g., lElmegreenl [l982l H9871: iKim et all [200l . In this 
subsection, we derive reduction factors for the conven- 
tional Jeans modes and the distortional modes in trun- 
cated polytropic disks. 

We first consider a pressureless, non- rotating disk (i.e., 

hi = re 2 , = 0) f° r simplicity. For given pi, equation 
has a formal solution 

2ttG 



pi(z')e~ 



\dz' 



(37) 



(e.g.. IKim fc Ostrikerll2007l) . Substituting equation 
into equation (f27|) . one obtains 



where 
F{k x 



£0X1 



p (z)p 1 (z')e~ 



(38) 



dzdz', (39) 



is the generalized reduction factor of self-gravity at 
wavenumber k x . It is apparent from equation (|39j) that 
T = 1 + 0(k x H) in the long wavelength limit, regardless 
of density distributions. 

Using equations ([2"5]) and ([2l)| , we decompose equation 
(|39|) into two parts as 



T(k x ) = wT } {k x ) + (1 - w)Tu(k x ), 



(40) 



where 



p Q (z)p 1 . l (z')e- k ^ z - z '\dzdz', (41) 
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Figure 5. Various gravity reduction factors as functions of k x H. 
The solid line plots Tj by ta king p o cx sech 2 (z/_ff) and pi from the 
eigensolutions of equations J20J— J23), while the dot-dashed line 
gives J-j by assuming pi oc po oc sech 2 (z/_ff). The reduction factor 
Jx> (eq. |45] ) of the distortional modes is plotted as a dashed line. 
For comparison, the commonly used T = (1 + kxH)" 1 is plotted 
as a dotted line. 

and 



•MM 



-k x a 



So 



Po(z) cosh(k x z) dz, (42) 



represent the reduction factors of self-gravity for the con- 
ventional Jeans modes and distortional modes, respec- 
tively. 

We integrate equation (|4"Tj) for an unbounded (A = 1), 
isothermal disk with po cx scch" (z / Hq) . The resulting 
J- j with pi.i cx po is plotted in Figure [5] as a dot-dashed 
line, while the case with p\^ numerically obtained from 
the full solution of equations ([2^|) -([2^ ]) is shown as a solid 
line. Note that both are closely approximated by (1 + 
k x H)~ x shown as a dotted line. For pressure-confined 
disks, we empirically found by varying A as well as the 
BCs that 



1 



and 



k x HA^I 21 
1 



for rigid BC, 



1 + KH ' 



for free BC, 



(43) 



(44) 



match the numerical dispersion relations of GI quite well 
(see below) Q The difference between equations (|43| and 

For an unbounded, exponential disk with po oc pi,; oc 
exp(— \z\/H), equation (14 1 1 yields 



T-CXp 

^3 



k x H 



l + k x H 2{l + k x H) 2 ' 



which is different from the often- used expression (eq. 1441 ). The lat- 
ter is in fact the exact force correction term at the z = plane (e.g., 
Kim & Ostriker 2007). while the former is the correcti on f actor av- 
eraged along the z-direction. Nevertheless, equation 1441 matches 
the reduction factor for sech 2 (z/H) disks better than ^"j Xp . 
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Figure 6. Dispersion relations of the fundamental modes (thick lines) and the first harmonics (thin lines) of the (a) acoustic and (b) 
acoustic-surfacc-gravity waves in disks with A = 0.8 and 7 = 1. Solid and dashed lin es p lot the numerical results with and without 
self-gravity, respectively, while dotted lines draw the approximate dispersion relation (eq. [46]) for the fundamental modes. 



(03]) is due to the fact that p\ * in equation (["4"Tj) is affected 
by the BCs. 

For small A, T is mostly from J"b (see eq. [40]). In this 
case, the distortional modes resulting from pi jS dominate, 
and it is reasonable to take po ~ constant and a w H in 
the evaluation of J-b • Equation (|42|) then becomes 

_ t „sinhfc T .ff 

(45) 

which is plotted as a dashed line in Figure[5] Incidentally, 
J-"d is not much different from (1 + k x H) , either. 

4.2. Non-rotating Disks 

We now explore GI of non-rotating, pressure-confined 
disks with the effect of thermal pressure included. Com- 
bining the results of Sections [3] and 14. 1[ we write an ap- 
proximate dispersion relation 

uj 2 = c 2 ctf kl - 2TrGZ k x (wFj + (I - w)F D ), (46) 

where c c ff is given by either equation (|30[) or equation 
depending on the BCs. 

4.2.1. Rigid Boundary 

First, we impose the rigid BCs with w = 1. Figure 
[5^ plots the dispersion relations of waves with (solid 
lines) and without (dashed lines) self-gravity in disks 
with 7 = 1 and A = 0.8 for the fundamental acoustic 
modes (thick) as well as the first harmonics (thin) ob- 
tained numerically by solving the perturbation equations 
(f!?0"|) - (P5]) . The approximate dispersion relation ([4*51) for 
the fundamental mode is plotted as dotted lines, which 
are in good agreement with the numerical results for 
k x tl <, 2. Clearly, self-gravity lowers the frequency and 
it is the fundamental mode that becomes unstable upon 
inclusion of self-gravity. GI under the rigid BCs is the 
conventional Jeans instability since surface distortion is 
absent. 

Figurc[7]plots as solid lines the numerical dispersion re- 
lations of pure Jeans modes for some selected values of A 



and 7, which agree quite well with equation (|46[) drawn as 
dotted lines: the relative differences between the true and 
approximate values of the maximum growth rates and 
critical wavenumbers are less than 5% for 0.8 < 7 < 2 
and < A < 1. For fixed So and c s q, the Jeans modes 
are more unstable for smaller A because the disk becomes 
geometrically thinner. In the limit of A — > 0, the disper- 
sion relations of GI in pressure-confined disks become 
identical to those of an infinitesimally-thin disk with the 
same So and c s o, shown as dashed lines in Figure [7K-c. 
For fixed A, disks with larger 7 are more unstable owing 
to a steeper temperature gradient as well as a smaller 
vertical scale height (H = Acl / (irGT, j) cx 7 -1 ). Fig- 
ure [Tji replots Figure "7^, by scaling oj and k x in terms 
of c e ff rather than the midplane sound speed c s q- Close 
agreement among the curves with different 7 evidences 
that the effective sound speed defined in equation (|30|) is 
really a good representative of the mean averaged sound 
speed of polytropic disks. 

4.2.2. Free Boundary 

Under the free BCs, perturbations experience two 
types of restoring forces: compressibility and surface 
gravity. The former and latter dominate when A ~ 1 
and A <C 1, respectively. For general A, these are mixed 
together to become acoustic-surface-gravity waves. Fig- 
ure [Bja plots the numerical (solid lines) and approximate 
(dotted lines; eq. [46]) dispersion relations for the fun- 
damental modes (thick) and the first harmonics (thin) 
in disks with A = 0.8 and 7=1. Non-sclf-gravitating 
waves are plotted as dashed lines, while solid lines draw 
the self-gravitating waves. The close agreement between 
the numerical and analytic results shows that the fun- 
damental modes under the free BCs are really acoustic- 
surface-gravity waves that become unstable when self- 
gravity is included. In what follow, we call the unsta- 
ble acoustic-surface-gravity modes the mixed GI. Figure 
[5] plots as solid lines the numerical dispersion relations 
of the acoustic-surfacc-gravity modes in non-rotating, 
pressure-confined disks with various values of A and 
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Figure 7. (a)— (c): Dispersion relations of the pure Jeans modes in non-rotating, pressure-confined, polytropic disks under the rigid BCs. 
The wavcnumbcr and growth rate are no rmalized by using the midplane sound speed. The solid lines are the numerical results, while the 
approximate dispersion relations (cq. [46]) arc plotted as dotted lines. The razor-thin dispersion relation (A — >■ 0) is shown as a long dashed 
line in each panel, (d): The A = 1 case plotted in (a) is redrawn using c e fj in the normalization of ui and k x . 



7. The wavenumber and growth rate are normalized 
by the effective half-thickness H and the free-fall time 
at the midplane, (27rGpoo) _1 , respectively. The cases 
with 7 = 1 recover the results of earlier studies (e.g., 
bimon 119651 lElmeereen fc Elmegreenlll978l lUsami et all 
1995t llwasaki et al.ll201lD . Shown also as dotted lines 
are the approximate dispersion relation in good agree- 
ment with the numerical results. Regardless of A, the 
maximum growth rate, |w max |, of unstable modes occurs 
at kH ~ 0.3-0.4. However, the physical nature of GI 
is markedly different depending on the degree of pres- 
sure confinement. When the external pressure is weak 
(Am 1), unstable modes are dominated by acoustic- type 
perturbations (w ~ 1), making the mixed GI under the 
free BCs similar to the pure Jeans mode occurring under 
the rigid BCs. When perturbations from surface distor- 
tion dominate (A <C 1), on the other hand, equation (|46[) 
is reduced to 



g 0a k x tanh^if) - 27rG£ fc :r x 



1 



-2k T H 



2k x H 
instability 



-, (47) 
(e.g., 



for the pure distortion al 
IGoldreich fc Lvnden-Belll[l965al) . 

Figure [9] compares the marginal wavenumber k x ^ c , the 
most unstable wavenumber k x , max , and the maximum 
growth rate oj max of the pure Jeans instability with those 
of the mixed GI of isothermal disks as functions of A. 
The mixed GI is identical to the pure Jeans mode when 
A = 1, but the former in general has a higher growth 



rate and occurs at a shorter wavelength than the lat- 
ter for arbitrary A. This is because the effective sound 
speed varies with A differently depending on the BCs. 
For A <C 1, c c ff ~ c s o for the pure Jeans modes (eq. 
[3~1] ) , while c cr T ~ y/goa,a oc A for the mixed GI (eq. [3S] ) . 
Since both types of GI have the characteristic wavenum- 
ber k XtC ~ 7rGSo/Ceff and growth rate \u>\ ~ k XtC /c c s, 
the critical wavenumber and maximum growth rate of 
the Jeans modes are smaller by a factor of A 2 and A, 
respectively, compared to those of the mixed GI. 

4.3. Rotating Disks and Stability Criteria 

Finally, we include the effect of disk rotation and study 
the stability condition to axisymmetric GI. In view of 
equation (|27[) . one may assume that the gas motions 
arising from thermal pressure and self-gravity are sep- 
arable from the epicycle motions, which is shown valid 
for k x H 1 in the preceding sections. The approximate 
dispersion relation of GI then becomes 



(48) 



which is a generalization of equation (fTJ) to pressure- 
confined, polytropic disks with vertical stratification. 
Figure [TU] compares equation (|4"B"|) shown as dotted lines 
with the full numerical results plotted as solid lines 
for disks with A = 0.7, 7 = 0.8, 1, 1.5, 2, and Kq = 
0.2(7rG£o/c s o) 2 j which arc in good agreement with each 
other, demonstrating again that equation (|48[) is accurate 
in describing GI of pressure-confined polytropic disks 
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Figure 8. Dispersion relations of the mixed GI in non-rotating, pressure-confined polytropic disks with 7 = 0.8,1,1.5,2 and A = 
0.4, 0.6, 0.8 under the free BCs. Solid lines draw the numerical results, while dotted lines plot the approximate dispersion relation (eq. 1461 ). 
In each panel, the pure distortional modes in the limit of A — > are shown as a long dashed line. 



provided c e fj and T are chosen appropriately. 

We generalize Toomre's stability parameter originally 
defined for a razor-thin disk to vertically-stratified, 
pressure-confined disks as 



/off 



(49) 



ttG£ 

Using c c ff/H = (2ajnGpoo) 1 ^ 2 for A = 1, one can show 
that the stability parameter Qm = ^ 2 /(47rGpoo) used by 
iMamatsashvili fc Ricd (I2010D is equal to \/Qcttl 
Since T < 1, the critical value of (Jeff is less than unity. 
Note that in the limit of strong pressure confinement, 
Qoff is reduced to 



nG^Q \/llGp 



for A ->• 0, 



(50) 



11(1 



applicable to the pure distortional modes. 

Figure [TT] plots the critical values Q c e,c for axisym- 
metric GI in the (7, A) plane. Figure Ilia is for pure 
Jeans instability under the rigid BCs, while the case 
of mixed GI with the free BCs is shown in Figure 
[Tib . As A — > 0, QeS, c always converges to 1, for pure 
Jeans modes, corresponding to a razor-thin disk, and 
to 0.756 for th e mixed GI, corresponding to a n incom- 
pressible disk (|Goldreich &: Lvnd cn-Bcll 19653). In the 
pure Jeans instability, finite disk thickness weakens self- 
gravity, decreasing <2 G ff,c with increasing A. For exam- 
le, an unbounded isothermal d isk has Q c ff, c = 0.676 
Goldreich fc Lvnden-Belll Il965al ). For the mixed GI, 
QeS,c ~ 0.68-0.75 largely insensitive to 7. 



5. SUMMARY AND DISCUSSION 

We investigate GI of rotating, pressure-confined, 
vertically-stratified gas disks to axisymmetric perturba- 
tions. As an initial equilibrium, we consider a self- 
gravitating, polytropic disk in hydrostatic equilibrium. 
The disk is truncated by a constant external pressure 
characterized by a dimensionless parameter A (eq. [TJ ) . 
To distinguish distortional modes of GI occurring in 
highly confined disks from the conventional Jeans modes, 
we adopt the rigid-surface BCs that allow only Jeans 
modes as well as the free-surface BCs under which two 
modes coexist. Our model disks do not include convec- 
tive motions. 

By deriving approximate dispersion relations and com- 
paring them with the numerical results from a set of 
the perturbation equations (|20|) -([23 j) . we find that GT 
unstable modes are in general even-symmetry fundamen- 
tal modes that have no vertical node. Under the rigid 
BCs, the fundamental modes are simply acoustic waves 
that propagate via thermal pressure only. Under the 
free BCs, on the other hand, surface-gravity waves aris- 
ing from surface distortion provide an additional restor- 
ing force, forming acoustic-surface-gravity waves. In the 
presence of self-gravity, the acoustic waves become unsta- 
ble to pure Jeans modes, while the surface-gravity waves 
become unstable to incompressible, distortional modes. 
Therefore, disks under the free BCs are unstable to the 
mixed GI in which the Jeans and distortional modes co- 
exist. When pressure confinement is weak (A ~ 1), the 
Jeans modes dominate, while the distortional modes be- 
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Figure 9. (a) Maximum (k x>m ^; dashed lines) and critical (k x ^ c ; solid lines) wavenumbers, and (b) squared maximum growth rate as 
functions of A for 7 = 1. The thick lines correspond to the mixed GI under the free BCs, while the thin lines are for the pure Jeans modes 
under the rigid BCs. 




Figure 10. Dispersion relations of rotating, polytropic disks with 
A = 0.7 and 7 = 0.8, 1, 1.5, 2 when k? q = 0.2 (7rGS /c s0 ) 2 . Solid 
lines plot the numerical results, whil e dotted lines draw the ap- 
proximate dispersion relation (eq. 1481 V 

come prevalent in strongly confined disks. The relative 
importance of the Jeans modes can be measured by the 
dimcnsionlcss parameter wo defined by equation (|35j) . 

While polytropic disks are vertically stratified both 
in density and temperature, we find that the effective 
sound speed c e ff defined either by equation (ptU|) for the 
rigid BCs or equation (|32|) for the free BCs represents 
the mean propagation speed of acoustic-surface-gravity 
waves quite well. Under the rigid BCs, c e g is given by 
the density- weighted, harmonic mean of the local sound 
speeds, while it has a contribution from the surface- 
gravity waves under the free BCs. 

We derive the reduction factors J-(k x ) responsible for 
the reduced gravity due to finite disk thickness for the 



pure Jeans modes (eqs. [43] and [44]) and for the dis- 
tortional modes (cq. |45|). The approximate dispersion 
relation (|48[) matches the numerical results very closely 
as long as the effective sound speed and the gravity re- 
duction factors are properly chosen. Disk rotation in- 
troduces inertial waves associated with epicyclic motion, 
but it is again the fundamental mode of acoustic-surfacc- 
gravity waves that are subject to GI. Using the effective 
sound speed, we define the stability parameter Q c g as in 
equation (|49|) that generalizes the Toomre's parameter 
into pressure-confined, polytropic disks. The pure Jeans 
modes have the critical values Q c ff,c = 0.67 - 1 with the 
larger value corresponding to a razor-thin disk, while the 
mixed GI has Q c g,c = 0.67 - 0.76 insensitive to 7. 

As mentioned in Introduction, the origin of the mixed 
GI of pressure-confined disks has been uncertain. Various 
mech anisms such as neutral modes (iLubow fc Pringle 
1993 ) , pinching force (lElmegreenl 119891 lUmekawa et al. 
1999[ iWunsch et al l |20JDD, and surface-gravity waves 
(jlwasaki et al.l l20l"TI ) have been proposed as the cause 
of GI. In this paper, we clearly show that the mixed GI 
occurring under the free BCs is a combination of pure 
Jeans and distortional modes, the latter of which is gen- 
erated by deformation of free surfaces. The existence of 
an additional restoring force due to the vertical gravity 
excites low-frequency fundamental modes, as opposed to 
high-frequency acoustic modes. With a reduced effective 
sound speed, the mixed GI has a larger growth rate and 
a smaller le ngth scale than the p ure Jeans modes. 

Recently, llwasaki et al.l (|2011| ) performed a linear sta- 
bility analysis of expanding isothermal shells around H II 
regions. They considered combinations of the rigid and 
free BCs as well as an asymmetric density profile caused 
by the bulk deceleration of a shell. They found that 
the numerically calculated dispersion relations are well 
described by an approximate relation lu 2 = c 2 s jk 2 — 
27rGI] fc x , where c 2 SI = [B x AttGT, H + (c s /2) 2 }. Al- 
though the dimcnsionlcss parameter B = 0.39 in their 
effective sound speed c e ff./ gives the best fit to their 
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Figure 11. Contours of the critical Q e g values of (a) the pure Jeans modes under the rigid BCs and (b) the mixed GI under the free 
BCs for pressure-confined, polytropic disks. As A — > 0, Q c ff lC approaches unity for the pure Jeans modes and 0.756 for the mixed GI, 
independent of 7. 



numerical growth rates, its physical nature is some- 
what ambiguous. Noting that H = A 2 c 2 / (ttGJ^q) , 
the above relation can be arranged as w 2 = c 2 s k x — 
2TiGYiQk x Fi{k x ) , similarly to equation (|4"6")l , with the 
effective sound speed of acoustic-surface-gravity waves 
c ff = 2A 2 /(1 + A 2 )c 2 in the long-wavelength limit 
(eq. [33]) and the gravity reduction factor modified to 
Fi{k x ) = 1 - k x H x [2B + (1 - 7A 2 )/(8A 2 (1 + A 2 ))}. 
For A > 0.4 (or J _» ex t/(7rGE ) Z 21/8) considered by 
llwasaki et ail (120111 ). 7i = 1 - (0.4-0.6)^7? in reason- 
able agreement with the reduction factors that we found 
(see Figure [5]) . This sugg ests that c Pr T.r in the approxi- 
mate dispersion relation of llwasaki et al.1 (|2011l ) contains 
a contribution from the reduced gravity due to finite shell 
thickness. 

While the mixed GI of pressure-confined disks has re- 
ceived a little attention compared to pure Jeans modes, 
it may be responsible for structure formation via frag- 
mentation of thin shells in the ISM. For example, 
iBovd fc Whitworthl (|2005l) claimed that GI of shock- 



compressed layers can explain the formation of Jupiter- 
mass (significantly smaller than t he Jeans mass) objects 
detec ted in you ng star cluster s dZapatero Osorio et ahl 
120021) . Recently. IWiinsch et al.l (120121) reported that the 
observed slopes of the mass spectrum of molecular clouds 
in the Carina Flare are indeed consistent with those ex- 
pected from GI of pressure-confined shells. Distortional 
modes in flattened systems can be thought of as aris- 
ing from the tendency to lower the gravitational poten- 
tial energy by trans forming into a spherical configuration 
(|Usami et al.l 119951 ). Thus, the mixed GI would occur 
not only in a planar geometry but also in a cylindri- 
cal gaseous column. In the presence of magnetic fields, 
distortional instability is su ppressed in the dir ection par- 
allel to the magnetic fields (|Nagai et al.lll998l) . naturally 
leading to cloud formation in filamentary shapes. 

We are grateful to the referee for a helpful report. This 
work was supported by the National Research Founda- 
tion of Korea (NRF) grant funded by the Korean gov- 
ernment (MEST), No. 2010-0000712. 



APPENDIX 

A. HYDROSTATIC EQUILIBRIUM 

In this appendix, we derive an analytic solution of equations ([7]) and ([5]) for the density distributions of self-gravitating 
polytropic disks. We begin by defining the dimcnsionless variables 

£ = z/H , (Al) 

e(0 - Wax)) 17 ", (A2) 

where n = 1/(7 — 1) is the polytropic index. Then, equations ([7]) and (|5J| result in 

Equation (|A3[) is the usual Lane-Emden equation in the planar geometry (e.g., iViala fc Horedtlll974D . which can be 
solved with the proper boundary conditions 8(0) = 1 and dO(0)/d £ — 0. 
Now, we combine equations (jHJ and (j9|) together with equations (|All) and (|A2j) to obtain 
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On the other hand, by integrating equation ((5J) over z and using equation (|A3[) . we express p in equation Q as 



where the second equality utilizes equation (|A4[) . The solution of equation (|A5|) that satisfies the boundary conditions 
is simply 

e = (i-/i 2 ) 1/(n+1) , (A6) 

or 

PoM=Poo(l-M 2 ) 1/7 , (A7) 

(see also. lHarrison fc Lakelll972h . 
To find po in terms of z, we substitute equations (|A1[) and (|A6[) into equation (|A4[) to obtain 



z /"^ dp 

A formal solution of equation (|A8I) is given by 



(A8) 



^= 2 F 1 (i,i;|;/, 2 ) M , (A9) 
where 2-P1 denotes the Gaussian hypergeometric function. Equation (| A8|) is intcgrable for a few particular values of 7: 
1. When 7 = 2/3, 

3/2 



z/ffo =/Vv / l-M 2 , and p = p 00 (l + (z/ff ) 2 ) , (A10) 



where Hq = K/(2ttGpq / q 3 ) (e.g., ILarsonllTMl . 

2. When 7=1, 

z/i?o = tanh -1 /i, and po = poosech 2 (z/7?o), (All) 
where = K/{2-KGp m ) (e.g., ILedomdH95lt IGoldreich fc Lvnden-Belll[T965ah . 

3. When 7 = 2, 

z/i?o = sin - /a, and po = Poo cos(z / Ho) , (A12) 
where Hi = K/{2txG) (e.g.. IGoldreich fc Lvnden^BeIl[T965al) . 

4. When 7 = 00, 

z/H = p 1 and p = Poo, ( Al 3) 

(e.g., lHarrison fc Lakdll972[ ). 

B. BOUNDARY CONDITIONS 

Equations (|2"0|) - (j2"3")) possess a reflection symmetry with respect to the z = plane, that is, invariant under the 
transformations z — > — z, £i z — >• — £i z , i/'i ~~ ^ — V4; ^1 — ^ V*i — ^ ^l- Thus, perturbations are in general a 
superposition of even-symmetry modes and odd-symmetry modes. Of these, odd-symmetric modes satisf ying = 
£■{?.(— z), h-\(z) = —h-\(—z), t^ -\ (z) = — -0i(— z), i/j[(z) = tJjU —z) are shown to be stable against GI (e.g. lSimonTll965l : 



lElmegreen fc Elmegreenlll978l : iMamatsashvihfe Rice 201Q). Therefore, we in this work consider only even-symmetry 



modes for which it is suffice to write only two conditions 

fi»(0) = ^i(O) = 0. (Bl) 

The other two conditions on hi(0) and ipi (0) are automatically satisfied from equations ff20|) — (f23 ]) . 

At the boundaries truncated by the external pressure, we consider two different types of BCs: the rigid and free 
conditions. In the rigid BC, the boundaries are assumed to be fixed with a vanishing Lagrangian displacement 

£l*L =o =0, for rigid BC, (B2) 



(e.g.. iVoitl [1981) . This condition may be appropriate to describe boundaries of a thin layer, as produce d by colliding 
clouds, since the shocked interfaces are stable against distortion (lUsami et~aII[T995t flwasaki et al.ll20ll . 

The free BC allows for a deformed Lagrangian surface on which the pressure is maintained constant to its initia l 
value pext (|Goldreich fc Lvnden-Belll965atlElmegreen fc Elmegreenlll978l : lLubow fc Pringlel[T99llIwasaki et al.ll2011| ). 



Since p cxt = p(a + £u\ ) = PoL = „+Pi L_ +(£u dp /dz) z = a to first order, one obtains 



(B3) 
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where go a = dipo/dz\ _ is the vertical gravity at the disk boundary (e.g.. iGoldreich fc Lvnden-Belllll965aD . 

In self-gravitating disks, another condition comes from applying the divergence theorem to the perturbed Poisson 
equation at the interface. Integrating equation (|18p over the volume of a thin shell placed at the interface gives 

= ^G P ^ lz \ z _ a , (B4) 

where the relation dipi/ dz\ z=a + = —k x ipi\ z=a is used since the gravitational potential outside the disk satisfies Laplace's 
equation. 

C. METHOD OF NUMERICAL INTEGRATION 

For a given set of disk parameters (7, A, k 2 ,), we wish to find a dispersion relation uj = u}(k x ) that satisfies the 
perturbations equations (|20j) - (|23[) subject to the two BCs and two even-symmetry constraints. Since the perturbation 
equations are linear in terms of the independent variables (£i z , hi, ipi, ^4), we are free to choose one of the variables 
arbitrarily. We thus fix hi\ z=a = 1 ( or £i z\ z=a = 1) and take trial values for ipi\ z=a and uj 2 , while obtaining £,i z \ z=a (or 
h\\ z=a ) and ip'i\ z =a from equations (|B2|) - (|B4|) depending on the adopted BCs. We then integrate equ ations (f2(3|) - ([23"|) 
from z = a to z = using the fourth-order R ung e-Kutta-Gill method (|Abramowitz fc Stegunlll972[ ). At the z = 
plane, we check the even-symmetry conditions (|B1[) : if these are not fulfilled, we go back to the disk boundary, change 
ipi\z=a and ui 2 slightly based on the Newton- Raphson technique, and repeat the integrations. Solutions with accuracy 
of 10 -4 , that is, with |£i z (0)/£i z (a/2)|, |'!/4(O)/'0i( a )l < 10~ 4 are obtained typically within less than three iterations. 

D. MODE CLASSIFICATION AND EFFECTS OF SELF-GRAVITY 

Here, we first derive local dispersion relations of waves in the WKB limit, and compare them with the numerical 
results. This will help us verify that acousti c modes become unstable in the presence of self-gravity. We then discuss 
our results in comparison with the results of iMamatsa shvili &; Ricel (|2010() who mistakenly claimed that gravity makes 
inertial modes unstable. 



KWl\z=a ~ —J— 
dz 




D.l. Local Analysis 

We seek for the solutions of equations ([20 ]) - ([23 )) for local waves with k x , k z ^> Hq 1 . Defining k 2 = —X\ l d 2 x\j dz 2 
with xi = h\ + if) 1, it is straightforward to derive 

LO* - L0 2 (pC 2 s k 2 + K 2 ) + K 2 Pc 2 k 2 z = 0, (Dl) 

wher e k 2 = k 2 + k 2 and (3 = 1 + ipi/hi = 1 — 4nG po/ \c 2 k 2 ) is the correction factor for self-gravity (e.g.. iChandrasekharl 
119611) . The corresponding eigenfunctions are 

= A ( ~f3ujc s k x }(u 2 - 4) \ e iC«-*+*--+*.«) , (D2) 
V ~(3c s k z /uj j 

where A is an arbitrary constant. 
For waves with |/3|c 2 fc 2 ^> Kq or k~^> k z , equation (|D1[) yields two approximate solutions 

tu 2 p = (3c 2 k 2 + kI (D3) 

= *l a lff z 2 - (D4) 
u I3c 2 k 2 + nl K ' 

It can be shown that the group and phase velocities of the waves associated with to 2 are parallel to each other, while 

those with u> 2 arc perpendicular. This demonstrates that the former is acoustic waves (p modes) boosted by epicycle 
motions, while the latter is inertial waves (r modes). Note that inertial waves require non- vanishing k z for propagation. 

It is apparent from equation (|D3|) that the local acoustic waves become gravitationally unstable provided self- 
gravity is sufficiently strong with (3 < ~K^/{c 2 k 2 ). On the other hand, the inertial waves usually have k z H 3> 1 so 
that ui 2 ~ Kpfc 2 /fc 2 regardless of f3, indicating that the character and frequencies of in ertial wave s are almost unchanged 
by self-gravity. Figure [12] demonstrates this by directly comparing equations (|D3j) and (|D4|) shown as dashed lines 
with the numerical dispersion relations plotted as dots obtained by integrating equations (|20 |) -(f23 | subject to the 
rigid BCs. An isothermal disk with A = 0.8 and k 2 ,/ (2irGpoo) = 0.2 is considered for the numerical calculations 
using the method described in Appendix [Cj Only the results for the fundamental acoustic modes (i.e., no node in the 
z-direction) and the lowest-order inertial modes with a single node arc shown. The inset in Figure 112b enlarges the 
regions with 0.7 < k x H < 0.8 an d 0.1 5 < uj 2 / (2irGpoo) < 0.21, where the acoustic and inertial modes cross each other. 
In plotting equations (|D3j) and (|D4j) . we utilize the eigenfunctions and their derivatives at z = to calculate (3 and 
k z . Note that agreement between the local and numerical dispersion relations is excellent for the non-self-gravitating 
case. Self-gravity makes j3 ^ — 0.4 for the acoustic modes, making them unstable for 0.12 < k x H < 0.63. On the 
other hand, the inertial modes have /3c 2 fc 2 //tg £ 38 and thus remain stable even in the presence of self-gravity. 
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Figure 12. Comparison of the local dispersion relations (dashed lines) with the numerical dispersion relations (dots) for acoustic and 
inertial modes in (a) non-sclf-gravitating and (b) self-gravitating cases. For the numerical results, an isothermal disk with A = 0.8 and 
Kq/ (2-kGpoo) = 0.2 is taken and the rigid BCs are adopted. In (a), the local dispersion relations are almost identical to the numerical 
results. In (b), the acoustic and inertial modes cross each other near k x H = 0.77, and the inset zooms in on the crossing regions. While 
self-gravity changes the acoustic mode dramatically, the inertial modes are almost intact. 

As in our present work, [Mamatsa shvili fc Ricd (|2010() also calculated dispersion relations of waves in vertically- 
stratified, polytropic, rotating disks. Except for the internal gravity waves arising from convective motions that are 
missing in our models, their numerical results are qualitatively similar to oursQ In interpreting their results, however, 
they argued that the influence of self-gravity is strongest on the inertial mode and that GI results from the first 
harmonics of the inertial mode. This erroneous finding seemingly results from the observation that the long- wavelength 
unstable branch of the dispersion curves are smoothly connected to the short-wavelength part of the inertial modes, 
which occurs at k x H ~ 0.77 in Figure [l"2"b. They claimed that the frequency ordering w 2 < Kq < u; 2 which is valid in 
the non-self-gravitating case holds true also for the the self- gravitating counterpart, so that t he lo west-f requ ency modes 
arc always the inertial modes regardless of k x . However, our local dispersion relations (eqs. (|D3[) and (|D4[1 ) show that 
this is correct only when /3 > 0, i.e., only for waves with relatively short wavelength^ For gravitationally unstable 
modes with (3 < — ^/(c^fc 2 ), w 2 < < a; 2 < Kq. In Figure [T2b. the lower-frequency mode at k x H < 0.77 is connected 
to the higher-frequency acoustic mode at k x H > 0.77, while the higher-frequency mode at k x H < 0.77 is connected to 
the lower-frequency inertial mode at k x H > 0.77. When the inertial and acoustic modes are mixed together by having 
similar frequencies and wavenumbers, it is quite ambiguous to tell which one is which near k x H = 0.77. 

D.2. Comparison of Eigenfunctions 

A clear way to distinguish the modes is to compare the corresponding eigenfunctions since they contain information on 
mode characteristics. Inertial modes, being incompressible in nature, should involve very weak density perturbations. 
This can be readily seen from equation (|D2[) since 



Ipi/pol _ - «o| Kok x 



\uix/c s \ \f3u>\c s k x j3c s kl 



«1, 



(D5) 



for inertial modes with k x /k z <C 1. On the other hand, acoustic modes with k x /k z 3> 1 have \pi/ po\/\ui x /c s \ as 
c s k x /\uj p \ ~ 0(1), implying that they rely on relatively large density perturbations. 

Figure [13] plots the vertical profiles of the numerically-calculated eigenfunctions Ui x > u iz, and px/po of the fun- 
damental acoustic modes (left panels) and the lowest-order inertial modes (right panels). The background state is 
the same as in Figure [T2l The most-unstable wavenumber k x H = 0.37, well away from k x H = 0.77 where the two 
modes are mixed, is chosen for all modes. The solid and dashed lines correspond to the cases with and without self- 
gravity, respectively. The normalization is such that u\ x /c s = 1 at z = 0. The non-self-gravitating acoustic mode has 
w 2 /(27rC?poo) = 0.41, k z = 0, u lx /c s = 1, u lz — 0, and pi/po = 0.57, independent of z, which are in fact exact solutions 
(see Appendix IE)) . When self-gravity is included, the acoustic mode becomes unstable with o; 2 /(27rGpoo) = —0.16; 
the corresponding eigenfunctions varies with z only sightly with k z H = 0.11, but they do not have any node. For 



2 Note that Mamatsashvili & Rice (2010) considered disks with 
7 > 1 and p e xt = 0, automatically truncated at z = z cu t (see cq. 
1131 ). With the imposed free BCs, they identified surface-gravity 
waves (/ mode), although the destabilizing effect of surface distor- 
tion is absent in their models. 



3 This is usually the case in stellar interiors since c 2 ~ irGpoR 2 
from the condition of hydrostatic equilibrium, with R denoting 
the stellar radius. Then, /S ~ 1 - 0.1(A/i?) 2 > for waves with 
wavelength A < R. 
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Figure 13. Comparisons o f the cigenfunctions for the acoustic modes (left panels) and the lowest-order inertial modes (right panels) with 
k x H = 0.37 shown in Figure 1X21 The solid and dashed lines correspond to the self-gravitating and non-sclf-gravitating modes, respectively. 
All cigenfunctions are normalized such that u\ x = 1 at z = 0. Note that the inertial modes involve very little density perturbations and 
are almost unaffected by self-gravity, while the acoustic modes have \ui z /ui x \ -C 1 since they propagate primarily along the horizontal 
direction. 



acoustic modes, |ui 2 /ui x | <C 1 for k z /k x <C 1, consistent with the prediction of the local analysis. The amplitude of the 
perturbed density for the acoustic modes is of order unity relative to the velocity perturbations. On the other hand, 
the inertial modes have a node at around z/H ~ 0.57-0.61, regardless of the presence of self-gravity. Gravity makes 
little change in u\ x and small changes in u\ z and p±. The corresponding eigcnfrcqucncics and vertical wavenumbers are 
uj 2 I (2ttG poo) = 0.1954,0.1949 and k z H = 2.424,2.421 in non-self-gravitating and self-gravitating cases, respectively. 
Note that \pi/ po\/\ui x /c s \ <C 1 for the inertial modes, just as expected. 
All of the above results suggest that the acoustic modes are strongly affected by self-gravity to be unstable when 



Gravitational Instability of Disks 17 
self-gravity is sufficiently strong, and that the incrtial modes arc not much influenced by self-gravity and remain stable. 

E. SOLUTIONS FOR WAVES IN NON-SELF-GRAVITATING DISKS 
When = Kq = 0, equations (f20|) and (f2l) become 

-T 1 = K h i - -4(7*1 - 9otiz), (El) 
dz bj- c z s 

(E2) 

az 



In equati on (|E1|) . the vertical gravity go should be considered as being external. Substituting equation (|E2[) into 
equation (|E1[) . we obtain 



Tz{^) + U- k *) p ° h ^°> (E3) 

where the equilibrium condition dpo/dz = — <7oPo/c 2 is used. Note that equation (|E3[) together with the rigid BCs and 
the symmetry constraint at the midplane is in the Sturm-Liouville form with eigenvalue 

fc 2 ^ 2 /c 2 -fc 2 , (E4) 

indicating that the fundamental modes have the lowest frequency (thus most susceptible to GI). 

For a time being, we limit to an isothermal disk with pa oc s ech ( z/Hq) for which analytic solutions of equatio 
can be derived. In terms of the dimensionless variable p, (cq. |Allj ). equation (|E3[) becomes 



d 2 v dv ( X 2 \ . . 

^-^- 2 ^ + { 2 -l^) V = ^ (E5) 

where v = (1 — p 2 ) 1 / 2 hi and A 2 = 1 — (k z H ) 2 . The above equation is of the associated Legendre type whose solution 
is given by 

v = C 1 P? + C 2 P^ x , (E6) 
where C\ and C2 are constants to be determined, and P X ±A is the Ferrers' associated Legendre function defined as 

±A/2 



if* = (A=FM)(jf^) , (E7) 
(IWhittaker fc Watsonl[l96l . 

The even-symmetry condition, dv/dp\^=o, at the midplane requires C\ = C2. After some manipulations, we find h\ 
in terms of z as 

h\ = sinh(z/i?o) sin(fc z z) — k z Ho cosh(z/ Hq) cos(k z z), (E8) 



where k z = y k z — H 2 and the proportionality constant C\ or C2 is omitted. The vertical wavenumber that satisfies 
the rigid BC, dh\j dz\ z=a = 0, is then given by 

^2 _ f 0, for fundamental mode, (E9) 

z - n ~~ \ n 2 n 2 /a 2 + 1/Hq, for harmonics with order n = 1,2, 3, • • ■ . ^ ' 



Now we consider general polytropic disks, and seek for fundamental- mode solutions of (|E3|) . While equation (|E3|) 
cannot be solved in a closed form for arbitrary 7, it has the simplest but important solution 

hi —> constant, (E10) 



in the long- wavelength limit regardless of the BCs. This corresponds to equation (|E8j) with k z = for isothermal disks. 
For waves with k x H <§; 1, gas motions are restricted most ly to the horizontal direction. In this case, the acceleration 
of the gas along the vertical direction (— w 2 £i z in eq. }E2| ) becomes relatively unimportant, resulting in dh\jdz w 0. 
As Figure [2] shows, h\ w constant is reasonably good for k x H ^ 1. 

Another limiting case is a strongly confined disk with 4<1, for which po(z) = cons tant and w 2 /c 2 <C k 2 , the latter 
of which shall be verified a posteriori. Then, the even-symmetry solution of equation (|E3|) is 

hi oc cosh(fc x z), (EH) 

which corresponds to equation (|E8[) in the limit of k z ik x and z/Hq <C 1. In this case, V ■ £1 = and the system is 
essentially incompressible. Note that equation (j34|) gives w 2 /c 2 < (g na H/c 2 )kl fc 2 , as expected. 
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Up to now we have ignored the effect of disk rotation, but it is a simple matter to show that the above results are 
valid also for rotating disks provided k z in equation (|E4[) is changed to 

(E12) 



which is identical to equation (|D1|) . This implies that the local solutions presented in Appendix |D] are exact for 
non-sclf-gravitating, isothermal disks. 
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